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THE  EFFECTS  OF  CONDUCTIVITY  ON 
HIGH-RESOLUTION  IMPULSE  RADAR  SOUNDING, 
ROSS  ICE  SHELF,  ANTARCTICA 

Rexford  M.  Morey  and  Austin  Kovacs 


INTRODUCTION 

Results  of  impulse  radar  sounding  studies  of  sea 
ice  have  shown  that  where  there  is  a  preferred  cur¬ 
rent  direction  under  the  ice  cover  the  crystal  struc¬ 
ture  of  the  ice  becomes  highly  ordered.  This  includes 
a  crystal  structure  with  a  preferred  horizontal  c-axis 
that  is  oriented  with  the  local  current  (Kovacs  and 
Morey  1978, 1979, 1980).  The  result  is  that  when 
electromagnetic  energy  is  radiated,  in  the  megahertz 
frequency  band,  from  a  dipole  antenna  in  which  the 
electric  field  (E-field)  is  oriented  perpendicular  to 
the  azimuth  of  the  c-axis  of  the  ice  crystals,  a  weak 
bottom  reflection  (or  none  at  all)  is  detected.  But 
when  the  E-field  is  oriented  parallel  with  the  c-axis 
there  is  a  strong  reflection  from  the  ice  bottom.  It 
has  also  been  found  that  the  frequency  dispersion 
of  anisotropic  sea  ice  varies  in  the  horizontal  plane. 
This  was  demonstrated  by  the  center  frequency  of 
the  reflected  signal  spectrum  from  the  ice  bottom, 
which  was  maximum  when  the  E-field  was  parallel 
to  the  preferred  c-axis  direction  and  minimum  when 
the  E-field  was  perpendicular  to  it. 

Coring  through  the  Ross  Ice  Shelf  at  Site  J-9 
(Fig.  1),  Zotikov  (1979)  found  the  shelf  to  be  416 
m  thick,  and  from  examining  the  ice  core  recovered 
noted  a  bottom  layer  composed  of  6  m  of  saline 
(sea)  ice.  He  also  discovered  that  the  sea  ice  struc¬ 
ture  consisted  of  long  ice  crystals  with  a  preferred 
horizontal  c-axis.  This  discovery  suggests  that  it 
should  be  possible,  using  the  CRREL  impulse  radar 
system,  to  detect  the  existence  of  the  sea  ice  layer 
at  Site  J-9.  And,  through  analysis  of  the  amplitude 
of  die  reflected  wavelet  from  the  sea  ice/water  inter¬ 
face  versus  antenna  E-fleld  azimuth  orientation,  it 


should  be  possible  to  determine  the  alignment  of  the 
current  at  the  ice/water  interface. 

We  have  used  the  CRREL  impulse  radar  success¬ 
fully  in  the  Arctic  to  profile  sea  ice  thickness  and  to 
investigate  sea  ice  anisotropy  arising  from  the  pre¬ 
ferred  horizontal  azimuth  alignment  of  the  c-axes  of 
the  ice  crystals  at  the  bottom  of  the  ice  cover  (Kovacs 
and  Morey  1978, 1979, 1980).  In  addition,  we  have 
used  the  system  in  the  Antarctic  to  investigate  the 
slope  and  lateral  extent  of  the  brine  infiltration  layer 
in  the  McMurdo  Ice  Shelf,  the  thickness  of  the  shelf, 
and  saline  ice  growth  under  the  floating  Koettlitz 
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Figure  1.  Map  of  the  Ross  Ice  Shelf  showing  Site  J-9 
and  the  McMurdo  Ice  Shelf  survey  locations. 


Glacier  ice  tongue  (Kovacs  and  Gow  197S,  1977, 
Gow  and  Kovacs  1979,  Kovacs  et  al.  1982).  Based 
upon  these  results  and  preliminary  radar  range  cal¬ 
culations,  we  concluded  that  the  radar  system  could 
also  detect  the  sea  ice  growth  layer  under  the  Ross 
Ice  Shelf  at  Site  J-9  and  the  preferred  c-axis  azimuth 
of  the  ice  crystals.  Since  the  alignment  of  the  crys¬ 
tal  c-axis  in  sea  ice  has  been  found  to  be  parallel  with 
the  current  at  the  ice  boundary  we  also  expected  to 
be  able  to  predict  the  alignment  of  the  current  at  the 
ice/water  interface.  If  the  radar  sounding  at  Site  J-9 
was  successful,  the  radar  data  were  to  be  correlated 
with  seawater  current  information  to  be  measured 
at  the  bottom  of  the  drill  hole  at  Site  J-9  by  Nor¬ 
wegian  investigators.  The  optimum  goal  was  to 
demonstrate  the  ability  to  profile  the  thickness  of 
the  sea  ice  growth  layer  over  a  larger  area  of  the  Ross 
Ice  Shelf  to  determine  where  sea  ice  is  growing  or 
ablating  and  where  the  glacial  ice  may  be  melting. 
This  information  is  important  in  understanding  the 
mass  balance  of  the  ice  shelf. 


PROFILING  SYSTEM 


frequency  content  of  the  signal  from  the  various 
interfaces.  The  time-of-flight  information  is  used  in 
calculating  ice  shelf  thickness.  The  reflected  ampli¬ 
tude  data  can  be  used  to  determine  the  azimuth  at 
which  the  maximum  reflection  occurred  at  the  sea 
ice/water  interface,  and  thus  to  find  the  under-ice 
current  alignment  (Kovacs  and  Morey  1978, 1979, 
1980). 


THEORETICAL  CONSIDERATIONS 

The  electromagnetic  properties  of  all  materials  are 
described  by  the  following  constitutive  parameters: 
effective  conductivity  ae,  effective  permittivity 
ee  =  er-e0,  and  effective  permeability  pe  = 
where  er  and  u,  are  the  relative  dielectric  constant 
and  relative  permeability,  e0  =  8.85  x  10~12  farads/ 
meter  (F/m),  and  (i0=4nx  10"7  henrys/meter  (H/m). 
Plane  wave  propagation  is  governed  by  the  complex 
propagation  constant  y,  which  may  be  determined 
from 

7  =  o  +/0  =  [~  w2pee  +  /ojpoJ%  (1) 


The  impulse  radar  system  used  has  been  described 
in  the  literature  previously  referenced.  In  short,  it 
consists  of  timing  electronics  which  clock  an  impulse 
generator  and  sampling  head.  The  impulse  source 
generates  60-  to  1 50- V  impulses  with  a  band  width 
of  nearly  200  MHz  at  a  selectable  repetition  rate  of 
12.5, 25,  or  50  kHz  and  a  pulse  width  of  3  to  10  ns 
(depending  upon  the  antenna  and  therefore  the  fre¬ 
quency  used).  The  impulse  is  fed  to  a  transmitting- 
receiving  antenna  which  radiates  an  electromagnetic 
(EM)  wavelet  into  the  underlying  material. 

Three  different  antennas,  with  center  frequencies 
of  20, 80  and  300  MHz,  were  used.  The  highest-fre- 
quency  antenna  was  for  high  resolution  profiling. 

The  lower  frequency  antennas  were  used  for  signal 
amplitude  reflection  measurements  from  the  sea  ice 
bottom  versus  antenna  E-fleld  orientation.  The  re¬ 
flected  radar  EM  wavelet,  for  example  from  the  ice 
shelf  bottom,  is  displayed  on  an  oscilloscope  and 
graphic  recorder.  The  signal  data  are  also  stored  on 
an  analog  tape  recorder  for  later  playback  and  data 
processing.  The  EM  signal  information  recorded  is 
the  travel  time  to  and  from  the  various  reflecting 
surfaces  and  the  voltage  amplitude  of  the  reflected 
EM  wavelet. 

In  the  laboratory  the  recorded  data  can  be  played 
back  and  analyzed  on  a  frequency  spectrum  analyzer. 
In  this  process  the  exact  time  of  flight  to,  for  example, 
the  glacial  ice/sea  ice  or  the  sea  ice/water  boundary 
can  be  obtained,  as  well  as  the  voltage  amplitude  and 


where 

(3  =  real  phase  constant 
a  =  attenuation  constant 
u  =  2ir/=  angular  frequency 
f  =  frequency 

As  the  radar  pulse  progresses  through  the  ice  with 
a  phase  velocity  vm  =  co/ft  it  suffers  an  exponential 
attenuation  with  an  attenuation  constant  a.  Thus, 
the  attenuation  A  in  decibels/meter  (dB/m)  is 

A -20  a  log,0(e) .  (2) 


The  pulse  is  partially  reflected  at  an  interface  between 
two  different  materials  (e.g.  ice  and  seawater),  the 
reflection  coefficient  p  being 


^2  ~  *>1 
T?j+7Jj 
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where  the  complex  characteristic  impedances  t }  of 
the  media  are  given  by 


Studies  of  the  electromagnetic  probing  of  snow 
and  ice  require  a  firm  understanding  of  the  complex 
dielectric  constant  of  the  materials  as  a  function  of 
frequency,  temperature,  and  impurity  content.  The 
relative  complex  dielectric  constant  for  ice  is  ade¬ 
quately  expressed  by  the  following  Debye  equation: 
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€roo  =  relative  high  frequency  dielectric  con¬ 
stant  =  3.15  ±0.05 

ers  =  relative  static  dielectric  constant  (tem¬ 
perature-dependent) 

t  =  relaxation  time  (temperature-dependent) 
°dc  =  zero  frequency  conductivity 

The  form  of  eq  5  is  presented  to  explicitly  evaluate 
the  effect  of  DC  conductivity  on  the  properties  of 
polar  ice,  since  for  pure,  solid  ice  there  is  no  DC 
conductivity.  The  real  effective  parameters  in  eq 
1,  ae  and  ee,  are  related  to  the  complex  constitutive 
parameters  by: 


a  =  o'  +  we"  (S/m) 


€e  =  e'  -  ~  =  e'(F/m)  0 

where  e  equals  the  real  part  of  the  dielectric  con¬ 


stant,  and  the  imaginary  part  of  the  conductivity 
o”  equals  zero  at  the  frequencies  of  interest.  The 
effective  conductivity  is  a  function  of  both  o',  which 
is  the  zero  frequency  conductivity,  and  e",  which 
is  the  dielectric  relaxation  loss  factor  due  to  the 
oscillation  of  dipolar  molecules  in  a  time-varying 
electric  field.  Thus  we  specifically  separate  the  two 
energy  dissipation  processes  in  ice,  namely  conduction 
and  relaxation.  The  DC  conductivity  of  ice  is  a 
measure  of  the  impurities  within  the  ice. 

The  relaxation  time,  r,  is  temperature-dependent 
and  is  also  dependent  upon  the  concentration  of 
impurities  in  the  ice  (Gross  et  al.  1978).  For  “pure” 
ice  r  is  well  approximated  by  (King  and  Smith  1981) 

r  =  r0  exp  (£77?  7)  (8) 

where 

Tq  =  5.3  x  10-»«  s 

E  =  activation  energy  *  13.25  kcal/mole 

R  -  gas  constant  *  1.9865  cal/K-mole 

Table  1  lists  some  representative  materials  and 
their  electromagnetic  parameters  at  a  frequency  of 
80  MHz.  The  conductivity  has  the  greatest  influence 
on  the  attenuation,  while  the  velocity  of  propaga¬ 
tion  in  the  material  um  is  mostly  a  function  of  the 
relative  dielectric  constant. 

The  maximum  range  of  an  ice-probing  radar  sys¬ 
tem  is  a  function  of  the  system  parameters,  the  tar¬ 
get  or  interface  parameters,  and  the  electromagnetic 


properties  of  the  material  being  probed.  The  radar 
equation  or  the  system  performance  factor  Q  is 


[Et‘Ek’Gt‘Gk  *1)**, 

T  -  g-- 

64  rr3/2!* 


(9) 


where 

Er  and  ER  =  transmitting  (T)  and  receiving 
(R)  antenna  efficiency 

GT  and  GR  =  transmitting  (T)  and  receiving 
(R)  antenna  gain 

um  =  velocity  of  propagation  in  the  medium 
g  =  backscatter  gain  of  target 
o  =  target  scattering  cross-sectional  area 
o  =  attenuation  coefficient  of  medium  (eq  1) 
L  =  range  to  target 
/=  frequency  of  signal 


or 

G=  10  log  [tt^J  00) 

where 

in  =  minimum  detectable  signal  power 
Pin  =  power  input  to  the  transmitting  antenna 

The  effective  conductivity  of  polar  ice  is  variable 
and  depends  on  temperature  and  impurity  content. 
Published  values  for  the  bulk  DC  conductivity  of 
polar  ice  range  between  7  x  10~6  S/m  and  5  x  10-s 
S/m  (Bentley  1977,  Robin  et  al.  1969).  Therefore, 
the  attenuation  in  polar  ice  may  be  between  6.44  x 
10-3  and  4.60  x  10~2  dB/m  at  80  MHz. 

Using  data  from  Table  1  and  eq  3  and  4  the  re¬ 
flection  coefficient  at  a  boundary  between  several 
different  materials  was  determined  (Table  2).  In 
Table  2  the  magnitude  of  the  reflected  signal  is 
seen  to  be  greater  at  a  glacial  ice/seawater  boundary 
than  at  a  glacial  ice/sea  ice  boundary.  A  glacial  ice/ 
seawater  boundary  would  therefore  be  detectable  at 
a  greater  distance  than  a  glacial  ice/sea  ice  boundary. 

Resistivity  measurements  made  by  Bentley  (1977) 
at  Site  J-9  indicated  that  the  bulk  DC  conductivity 
was  1 .4  x  10"*  S/m.  Prior  to  our  field  trip  this  in¬ 
formation  was  used  to  estimate  the  maximum  depth 
of  penetration  (using  eq  9  and  10  to  estimate  maxi¬ 
mum  L)  for  the  radar  antenna  systems  we  used  (to 
be  discussed  later).  The  calculated  maximum  depth 
of  penetration  was  greater  than  600  m  with  and  with¬ 
out  sea  ice  on  the  bottom  at  sounding  frequencies 
of  20  and  80  MHz.  Thus  we  expected  to  see  the 
bottom  of  the  ice  shelf  at  Site  J-9. 
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Table  1.  Approximate  electromagnetic  parameters  at .80  MHz. 


Material 

ae 

(S/m) 

t*r 

A 

(dB/m) 

um 

(m/ns) 

n 

(ohms) 

Air 

0 

1 

1 

0 

0.30 

377 

Polar  ice 

10-5 

1 

3.15 

9.2  X  10"3 

0.168 

21 1.7+/ 0.075 

Sea  ice 

lb-1 

1 

6 

42.8 

0.078 

62  +  /  48 

Seawater* 

4 

1 

81 

295 

0.014 

9.2+/  8.5 

Fresh  water* 

10"3 

1 

81 

0.18 

0.033 

41.9+/ 0.058 

*  At  20°C. 


Table  2.  Reflection  at  a  smooth  boundary. 


Reflector  boundary 

Refl.  coeff.  (p) 

Refl  magnitude 

Phase  angle  (p) 

Air-ice 

-0.28+/ 0.0002 

28% 

180° 

Shelf  ice-seawater 

-0.91  +/ 0.073 

92% 

175° 

Shelf  ice-sea  ice 

-0.50  +  / 0.26 

56% 

152° 

FIELD  PROGRAM 

In  preparation  for  the  Site  J-9  field  program  sev¬ 
eral  changes  and  modifications  were  made  to  the 
CRREL  impulse  radar  system  in  order  to  improve 
its  performance  and  thus  increase  the  maximum 
radar  range.  The  antennas,  each  operating  at  a  dif¬ 
ferent  center  frequency  (20, 80  and  300  MHz),  were 
tested.  A  preamplifier  was  placed  at  the  front  ends 
of  the  receivers,  providing  an  additional  10  to  1 5  dB 
of  gain.  Generally,  this  radar  system,  like  most  other 
radar  subsurface  sounding  systems,  uses  one  antenna 


for  transmitting  and  receiving  simultaneously.  In 
our  study  two  sets  of  matched  dipole  antennas  op¬ 
erating  at  either  80  or  300  MHz  were  used,  both 
transmitting  and  receiving  simultaneously  (Fig.  2). 
Theoretically  this  would  provide  a  6-dB  gain  in  sys¬ 
tem  performance.  Experimental  results  indicated 
about  a  4.5-dB  improvement. 

The  20-MHz  system  consisted  of  special  liquid- 
filled  antennas,  one  for  transmitting  and  one  for 
receiving  (Fig.  3).  Each  antenna  consisted  of  a 
plastic  tube  10  cm  in  diameter  and  about  3  m  long. 
A  uniquely  configured  copper  dipole  element  was 
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Figure  3.  Twenty-MHz  transmitting  and  receiving  antennas  at  Site  J-9. 
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Figure  4.  Maximum  ice/water  interface  depth  detectable  as  a  function  of 
the  effective  (bulk)  conductivity  of  an  ice  shelf  having  a  smooth  bottom. 

Table  3.  Radar  parameters  for  calculating  maximum  range. 

Radar 
parameter 

'in  (peak  watts) 

^mln  <watts> 

Q 


20  MHz 


80  MHz 


300  MHz 


2.5  X  103 
2.5  X  lO"8 
-UOd 
3d 
2  d 


50 

5  X  10-** 


12 

1.2  x  10"11 


supported  along  the  center  axis  of  the  plastic  tube, 
which  was  divided  into  several  chambers.  A  low- 
loss  liquid  Tilled  one  chamber  and  a  high-loss  liquid 
the  other.  The  purpose  of  the  dielectric  liquid  was 
to  make  the  antenna  electrically  longer  than  its 
physical  length  and  thus  radiate  lower  frequencies. 
The  high-loss  dielectric  liquid  tended  to  damp  out 
pulse  ringing  on  the  antenna. 

Table  3  lists  the  radar  parameters  used  for  calcu¬ 
lating  the  maximum  radar  range.  The  80-  and  300- 
MHz  system  values  in  Table  3  are  based  upon  the 
dual  transceiver  antenna  configuration.  The  maxi¬ 
mum  radar  range  as  a  function  of  conductivity  for 
the  20-  and  80-MHz  antennas  is  plotted  in  Figure  4. 
Two  models  are  considered,  one  without  any  sea  ice 
on  the  ice  shelf  bottom  and  one  with  6  meters  of 
sea  ice.  The  maximum  radar  range  is  shown  to  be 
very  sensitive  to  the  effective  conductivity  of  the 
shelf  ice. 

Before  leaving  for  Site  J-9  we  tested  the  radar 
system  on  the  McMurdo  Ice  Shelf.  We  were  able  to 
penetrate  about  2S0  m  of  the  shelf  and  receive  fairly 
strong  ice/water  reflections  at  all  three  radar  sounding 
frequencies.  However,  the  weakest  bottom  return 
occurred  at  300  MHz.  Measurements  were  also  made 
on  the  McMurdo  Ice  Shelf  using  a  single  80-MHz 
antenna  suspended  from  a  helicopter  flown  at  an 
altitude  of  10-1 5  m.  Shelf  thicknesses  of  up  to  27S 
m  were  continuously  profiled  (Fig.  5).  This  was  the 
maximum  shelf  thickness  encountered  on  the  flight. 

The  ice  shelf  thickness  at  Site  J-9  is  416  m  (Zoti- 
kov  1979).  Radar  profile  traverses  extending  up  to 
2  km  from  Site  J-9  were  made  with  both  the  20-MHz 
and  the  dual  80-MHz  antennas.  Neither  antenna  sys¬ 
tem  appeared  to  detect  the  glacial  ice/sea  ice  bound¬ 
ary  or  the  sea  ice/seawater  boundary. 


DISCUSSION 

Since  our  earlier  estimates  indicated  that  we 
should  have  “seen”  the  bottom  of  the  ice  shelf  at 


Site  J-9,  we  reevaluated  the  radar  range  equation 
(eq  9)  and  the  parameters  that  contributed  to  it, 
and  also  used  more  recent  information  (Jezek  1980) 
on  the  attenuation  characteristics  of  the  ice  shelf  at 
Site  J-9. 

Using  the  temperature-depth  profile  at  Site  J-9, 
Jezek  (1980)  determined  the  two-way  absorption  loss 
to  be  1 5  dB  ±  1  at  35  MHz.  Equation  5  was  used  to 
calculate  e"  as  a  function  of  temperature  for  “pure” 
ice  (oDC  =  2.1  x  10"8  S/m)  at  35  MHz,  where  the 
values  for  ers  and  t  were  taken  from  King  and  Smith 
(1981).  The  effective  (bulk)  conductivity  oe  was 
then  calculated  using  eq  6  for  “pure”  ice  and  for  the 
ice  shelf  at  Site  J-9  with  a  bulk  DC  conductivity  of 
1.4  x  10"*  S/m  as  found  by  Bentley  (1977).  Table  4 
lists  the  ers  and  r  vs  temperature  values  used  in  the 
analysis  along  with  the  calculated  oe  results.  Figure 
6  is  a  plot  of  ae  at  35  MHz  as  a  function  of  tempera¬ 
ture.  For  the  416-m-thick  ice  shelf  at  Site  J-9,  a  two- 
way  attenuation  of  1 5  dB  as  determined  by  Jezek 
(1980)  represents  an  effective  (bulk)  conductivity  of 
1.9  x  10"*  S/m  (assuming  no  impurities  in  the  ice 
shelf)-  And,  for  a  DC  conductivity  of  1.4  x  10"* 

S/m  as  found  by  Bentley  (1977)  at  Site  J-9,  the  effec¬ 
tive  (bulk)  conductivity  of  the  ice  shelf  is  estimated 
to  be  3.3  x  10-*  S/m  at  35  MHz  (Fig.  6).  Note  that 
both  conductivity  and  relaxation  time  are  concentra¬ 
tion-dependent  (Gross  et  al.  1978).  Therefore,  the 
difference  in  the  two  oe  values  (1.9  x  10"*  vs  3.3  x 
10"*  S/m)  given  above  may  be  due  to  impurity  effects. 
Since  we  do  not  know  the  type  of  impurities  or  their 
concentrations  at  Site  J-9,  we  cannot  include  this 
effect  on  the  relaxation  time.  However,  impurities 
within  the  icc  tend  to  decrease  r  (increasing  e"),  and 
thus  would  increase  oe  beyond  the  value  derived  from 
Jezek’s  data. 

Figure  7  is  a  plot  of  effective  conductivity  as  a 
function  of  frequency  for  two  different  temperatures. 
Above  100  kHz  the  effective  conductivity  is  frequency- 
independent  but  is  dependent  on  temperature  and 
impurity  content.  If  one  assumes  that  the  average 
temperature  of  an  ice  sheet  such  as  the  Ross  Ice 


Table  4.  Static  dielectric  constant,  relaxation  time,  and  imaginary  part  of  relative  die¬ 
lectric  constant  of  pure  ice  (o^  *  2.1  x  10~8  S/m)  at  3$  MHz. 

T  T  ae  (S/m)  oe  (S/m) 

( Q  tft  (,) _ «;  (opc  ■ 21  *  io~s)  <°dc  - 14  * 

0  91.1  2.05  xl0-s  1.95  x  10"J  3.8x10"*  5.20x10-* 

-5  93.4  3,25  x10-*  1.26x10"*  2.46x10"*  3.85x10** 

-10  95.0  5.13  x10**  8.15  x  10-*  1.59x  10"*  2.99x10"* 

-15  96.1  8.70x  10-*  4.87  xlO"3  9.5  x  10*6  2.35  x  10"* 

-20  97.2  1,49  xlO"4  2.89  x  1  O'*  5.6  xlO"6  1.96x10"* 

-25  98.2  2.55  xlO"4  1.71x10"*  3.35  X 10"6  1.73x10"* 


0  -5  -10  -15  -20  -25 


Temperature  (*C) 

Figure  6.  Plot  of  effective  (bulk)  conductivity  at  35 
MHz  as  a  f  unction  of  temperature. 


Shelf  is  constant  over  large  areas,  then  the  effective 
conductivity  is  essentially  a  function  of  impurity 
content. 

From  the  absorption  loss  information  of  Jezek 
(1980),  we  have  now  made  an  estimate  of  the  appar¬ 
ent  high  frequency  (>  100  kHz)  effective  (bulk)  con¬ 
ductivity  of  3.3  x  10"s  S/m  for  the  ice  shelf  at  Site 
J-9.  This  value  was  determined  by  correcting  the 
bulk  DC  conductivity  of  Bentley  (1977)  for  tempera¬ 
ture  and  frequency  effects  as  depicted  in  Figure  7. 
Using  oe  -  3.3  x  10-5  S/m  we  can  now  recalculate 
the  maximum  sounding  depth  of  our  20-MHz  radar 
system  at  Site  J-9  to  be  about  425  m  with  no  sea  ice 


layer  on  the  ice  shelf  bottom  and  about  405  m  with 
6  m  of  sea  ice  on  the  bottom  (Fig.  4).  The  above 
analysis  indicates  that  our  radar  equipment  was  just 
at  the  limit  of  being  able  to  detect  the  shelf  bottom 
at  Site  J-9.  The  maximum  calculated  depth  of  pene¬ 
tration  for  the  300-MHz  system  is  about  260  m  for 
an  ice  shelf  having  a  at  =  3.3  x  10~5  S/m  and  a  smooth 
ice  bottom  with  no  sea  ice.  Thus  the  weak  reflection 
obtained  from  the  ice/water  interface  on  the  McMurdo 
Ice  Shelf  at  a  depth  of  about  250  m,  as  previously 
discussed,  was  apparently  due  to  the  fact  that  this 
system  was  approaching  its  maximum  sounding  depth. 

At  Site  J-9  we  did  observe  many  internal  layers 
and  parabolic  reflections  believed  to  be  signatures  of 
buried  crevasses  (Fig.  8).  The  crevasses  are  believed 
to  have  developed  upstream  of  Site  J-9  and  to  have 
become  buried  by  subsequent  snow  accumulations 
during  the  years  they  were  moving  toward  the  site. 

The  tops  of  these  relic  crevasses  are  calculated  to  be 
about  46  m  below  the  snow  surface.  Figure  9  shows 
concave  undulations  in  the  internal  layers  above  and 
between  relic  crevasses.  Clough  (1975)  and  Jezek 
(1980)  reported  on  bottom  crevasses  in  the  Ross  Ice 
Shelf  and  at  Site  J-9  using  radio-echo  sounding  tech¬ 
niques.  Because  of  the  length  of  the  pulse  used  in 
their  radar  systems,  near-surface  details  were  lost. 
Surface  topography  at  Site  J-9  does  not  indicate  the 
presence  of  near-surface  relic  crevasses  or  bottom 
crevasses.  However,  Jezek  (1980)  suggested  that  as 
some  bottom  crevasses  form,  a  surface  depression 
develops  between  them.  He  speculated  that  the  shelf 
surface  between  the  water-filled  bottom  crevasses 
would  settle  by  about  13  m,  forming  a  concave  de¬ 
pression.  Our  radar  profiles  at  Site  J-9  support  this 
concept,  and  we  propose  that  surface  crevasses  may 
have  formed  as  a  result  of  bottom  crevassing,  as  de¬ 
picted  in  Figure  10. 

As  pointed  out  earlier,  the  maximum  detectable 
bottom  depth  for  a  radar  sounding  system  is  in  part 


urfoce 


re  8.  Graphic  record  of  radar  profile  over  two  of  the  many  “side  by  side  ”  buried  crevasses  observed  near  Sit e  3-9.  Time  between 
:ontaI  scale  lines  is  100  ns. 


turn  not  shown.  Time  between  horizontal  scale  lines  is  100  ns. 


Figure  JO.  Schematic  of  proposed  method  in  which 
some  surface  crevasses  are  formed  as  a  result  of  bot¬ 
tom  crevassing. 


a  function  of  the  effective  conductivity  of  the  ice 
(Fig.  4).  The  conductivity  of  the  ice  is  a  function 
of  its  temperature  and  impurity  content.  Many 
radar  sounding  profiles  obtained  from  aircraft  flights 
over  the  Ross  Ice  Shelf  have  sections  of  data  showing 
an  absence  of  bottom  echo  (Jezek  1980).  Neal 
(1982)  evaluated  small-scale  roughness  on  the  bottom 
of  the  Ross  Ice  Shelf  and  its  effect  on  radio  echo 
fading  patterns.  The  lack  of  ice  shelf  bottom  return 
is  also  attributed  to  such  things  as  brine  infiltration, 
crystal  fabric  alignment,  and  debris  within  the  ice. 

We  suggest  that  bottom  echo  fading  and  extinction 
in  radio  echo-sounding  may  also  be  due  to  changes 
in  the  bulk  DC  conductivity  of  the  ice  shelf.  Bentley 
(1977)  pointed  out  that  measurements  near  Roose¬ 
velt  Island  and  at  Byrd  Station  indicate  a  DC  con¬ 
ductivity  one-half  that  at  Site  J-9,  or  about  7  x  10"6 
S/m.  This  represents  an  effective  (bulk)  conductivity 
of  about  1.6  x  10~s  S/m  at  our  radar  sounding  fre¬ 
quencies.  Referring  to  Figure  4,  at  the  latter  con¬ 
ductivity  the  maximum  depth  of  penetration  of  our 
radar  system  operating  at  20  MHz  would  be  about 
720  m  instead  of  about  425  m  at  Site  J-9  where  the 
apparent  effective  (bulk)  conductivity  is  3.3  x  10-* 
S/m.  Different  radar  systems  will,  of  course,  have 
different  performance  characteristics  and,  therefore, 
the  curves  in  Figure  4  will  be  shifted  vertically  up  or 
down  depending  on  the  system  used.  However,  rel¬ 
atively  small  changes  in  bulk  conductivity  will  change 
the  amplitude  of  shelf  bottom  echos. 

The  question  then  arises  as  to  what  mechanism 
causes  variations  in  conductivity  in  the  ice  at  dif¬ 
ferent  locations.  One  explanation  is  that  some  dif¬ 
ference  in  the  densiflcation  or  strain  history  of  the 
ice  may  causa  variations  in  conductivity,  since  it  is 
assumed  that  the  snow  that  fell  over  large  areas  con¬ 
tained  about  the  same  impurities.  However,  there 
is  evidence  that  accumulation  rates  and  micropar¬ 
ticle  distribution  vary  considerably  over  regions  of 
the  Antarctic  (Thompson  and  Mosley-Thompson 


1982).  Jezek  (1980)  pointed  out  that  the  amplitude 
of  bottom  reflections  seen  in  airborne  radar  sounding 
data  taken  on  the  Ross  Ice  Shelf  changes  from  strong 
to  very  poor.  For  example,  strong  bottom  returns 
are  seen  for  the  ice  emerging  from  East  Antarctic 
outlet  glaciers  feeding  the  Ross  Ice  Shelf,  whereas 
weak  radar  reflections  are  seen  for  the  ice  shelf  bot¬ 
tom  in  areas  between  these  outlet  glaciers.  In  addi¬ 
tion,  the  radar  bottom  reflections  diminish  in  inten¬ 
sity  as  the  ice  moves  downstream  toward  the  Ross 
Sea.  A  possible  implication  is  that  the  conductivity 
of  the  ice  is  less  in  the  outlet  glaciers  and  that  local 
storms  deposit  more  sea  salts  near  the  edge  of  the 
Ross  Ice  Shelf,  increasing  the  shelf  conductivity 
locally.  Jezek  also  mentions  the  existence  of  a  few 
strong  multiple  bottom  reflections  near  Roosevelt 
Island -stronger  than  those  near  Site  J-9.  The  re¬ 
ported  bulk  DC  conductivity  at  Site  J-9  was  about 
twice  that  near  Roosevelt  Island.  Thus,  the  difference 
in  radio-echo  strength  may  be  due  to  local  and  re¬ 
gional  changes  in  the  bulk  DC  conductivity  of  the 
ice  due  to  variations  in  impurity  content. 
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